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ABSTRACT 

The 6 Li abundance observed in metal poor halo stars appears to exhibit a plateau as a function 
of metallicity similar to that for 7 Li, suggesting a big bang origin. However, the inferred primordial 
abundance of 6 Li is ^1000 times larger than that predicted by standard big bang nucleosynthesis 
for the baryon-to-photon ratio inferred from the WMAP data. Also, the inferred 7 Li primordial 
abundance is 3 times smaller than the big bang prediction. We here describe in detail a possible 
simultaneous solution to both the problems of underproduction of 6 Li and overproduction of 7 Li in 
big bang nucleosynthesis. This solution involves a hypothetical massive, negatively-charged leptonic 
particle that would bind to the light nuclei produced in big bang nucleosynthesis, but would decay long 
before it could be detected. We consider only the Y-nuclear reactions and assume that the effect of 
decay products is negligible, as would be the case if lifetime were large or the mass difference between 
the charged particle and its daughter were small. An interesting feature of this paradigm is that, 
because the particle remains bound to the existing nuclei after the cessation of the usual big bang 
nuclear reactions, a second longer epoch of nucleosynthesis can occur among Y-nuclei which have 
reduced Coulomb barriers. The existence of the X~ particles thus extends big bang nucleosynthesis 
from the first few minutes to roughly five hours. We confirm that reactions in which the hypothetical 
particle is transferred can occur that greatly enhance the production of 6 Li while depleting 7 Li. We 
also identify a new reaction that destroys large amounts of 7 Be, and hence reduces the ultimate 7 Li 
abundance. Thus, big-bang nucleosynthesis in the presence of these hypothetical particles, together 
with or without an event of stellar processing, can simultaneously solve the two Li abundance problems. 
Subject headings: cosmological parameters — dark matter — early universe — elementary particles 
— nuclear reactions, nucleosynthesis, abundances — stars: Population II 



1. INTRODUCTION 

A long-standing effort in astrophysics has involved 
searches for signatures of unstable particles that might 
have existed in the early Universe, but have long since 
become extinct. The precision that exists in the abun- 
dances of the nuclides produced in big bang nucleosyn- 
thesis (BBN) suggests that primordial nucleosynthesis is 
a good place to look for such signatures. In that context, 
it is of considerable i nterest that recent observations (e.g. 
lAsplund et ail 120061 ) of metal poor halo stars (MPHSs) 
indicate that the primordial abundances of both 6 Li and 
'Li may not be in agreement with the predictions of 
standard BBN. Specifically, the 6 Li appears to have an 
abundance plateau similar to that for 7 Li in very low 
metallicity stars. This suggests a primordial abundance. 
However the abundance value is roughly a factor of 1000 
larger than predicted in standard BBN. A simi lar, though 
much less severe, problem exist s for 7 Li (|Ryan et al.l 
[20051: iMelendez fc Ramirez! [200l lAsplund et all l2006f ): 
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the BBN value that would be consistent with the baryon- 
to-photon ratio inferred from analysis of the Wilkinson 
Microwave Anisotropy Probe (WMAP)-power spectrum 
suggests an abundance that is roughly a factor of three 
higher than the observed value. 

In response to this intriguing 6 Li result, a number 
of solutions have been suggested. Some relate the Li 
anomalies to the possible exis tence of unstable heavy pa r- 
ticles in the early Universe (Dimopoulos et al. 1 988alfbl . 
| 1989HJedamzik|2000l:lCvburt eT al. 2003: Jeda mzild[200l 
Kawasaki et al.l 120051 : iKusakabe et alj I2006D which lead 
to 6 Li production via non-thermal reactions induced by 
the particle decay. Some have suggested a possible epoch 
of enhan ced cosmic ray nucle o synthesis via the a + a re- 
actio ns (|Rollinde et alj 120061 : iPrantzosI [2001 iKusakabd 
[2003). 



Of interest 
sugge s tion 



to the present work, ho wever, is the 



recent works 



(ICvburt et all 
""20071: 



20071 : 



several 

2006; Kaplinghat & Raia rama nl 12001 [Pospelovl 
Kohri fc Takavamal [20071: lHamaguchi et al.l 
Kusakabe et alj 120071 : iJedamzikl l2008allbT ) that heavy 
negatively charged unstable particles could modify BBN 
and lead to 6 Li production. However, this mechanism 
would operate in a rather different way from the previous 
suggestions. In this paradigm, the heavy particles (here 
denoted as X~ ) would bind to the nuclei synthesized 
during BBN to produce exotic nuclei (hereafter denoted 
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as X-nuclei). These massive X~ particles would be 
bound in orbits that would have radii comparable to 
those of the nuclei to which they were attached. Their 
presence, however, would reduce the Coulomb barrier 
seen by incident charged particles, thereby enhancing 
the thermonuclear reaction rates and allowing a longer 
duration of BBN in an extended epoch of nucleosynthe- 
sis. In the present work we expand upon our previous 
study (jKusakabe et al.l l2007h and describe details of a 
revised analysis of X~ particle effects on BBN. The 
focus of the present work is to better quantify the effects 
of X~ particles on BBN. It is possible, however, that 
both X~ effects and the X-particle decay contribute to 
the abundances. In a future work we will consider this, 
but for now we only address the effects of X~ particles 
and show that this alone can explain the observed 
lithium abundances. 

In standard BBN, 6 Li production is suppressed. It 
is synthesized primarily via the 4 He(d,7) 6 Li reaction, 
which has a very small cross section. However. IPospeiovl 
(|2007f ) suggested that a large enhancement of the 6 Li 
abundance could result from an X~ bound to 4 He (de- 
noted 4 Hex). This allows for the X~ transfer reaction 
of 4 Hex(d,X~) e Li, which can enhance the production of 
6 Li by as much as seven orders of magnitude. Subse- 
quently, iHamaguchi et al.1 (|2007h carried out a theoret- 
ical calculation of the cross section for this X~ trans- 
fer reaction in a quantum three-body model. Their 
v alue was about a n order of magnitude smaller than that 
of IPospeiovl (|2007f ). This difference can be traced to their 
exact treatment of the quantum tunneling in the fusion 
process and the use of a better nucle ar potential. 

In other work. ICvburt et al.l (|2006l ) identified the X" 
as the supersymmetric counterpart of a tau, i.e., a stau., 
thereby providing some plausibility for the existence of 
such particles. They also considered the X~ transfer re- 
actions for 6 Li, 7 Li, and 7 Be produ c tion a s Pospelov sug- 
gested. iKaplingha t fc Raiaramanl (|200(i) observed that 
the decay of an X~ when bound to 4 He occasionally 
knocks out a proton or neutron to produce 3 He or 3 H, 
thereby enhancing their abundances in BBN. Those nu- 
clides could then interact with other primordial He nu- 
clei at higher energies than those normally associated 
with the B BN production of 6 Li . iKohri fc Takavainal 
(|2007t) and iKusakabe et all (|2007l ) studied the recombi- 
nation of nuclei with X~ particles, and suggested the 
possibility that BBN with charged massive particles also 
destroys 7 Li and so solves the 7 Li problem as well. In 
addition, a resonant reaction r Bex(p,j) 8 &x has recently 
been proposed (Bird et al. 2007) that further destroys 
7 Bex through an atomic excited state of 8 Bx- 

In the present work we have estimated ne w reaction 
rates f or the X~ transfer reactions suggested bv lPospelovl 
(2007) as well as others that could occur. We have 
re analyzed the X~ t ransfer reaction rates estimated 
in ICvburt et al.1 (|2006l ) and found them to be negligible 
based upon the reaction dynamics discussed below. 

We have adopted a simple model to estimate the bind- 
ing energy due to the bound X~ particles for each nu- 
cleus involved in BBN. The reaction rates are corrected 
for the modified nuclear charges and the effective mass 
resulting from the binding of the X~ particle. The re- 
action rates, however, are not particularly sensitive to 
the mass of the X~ particle. Therefore, we have simply 



assumed that they are much heavier than the nucleon 
mass. 

In this paper we study the potential consequences of 
X~ particles on BBN, including recombination of the 
X~ particles by the existing normal nuclei. Several new 
reaction processes have been included that can produce 
or destroy 6 Li and 7 Li in the early universe. We uti- 
lize a fully dynamical description which is solved numer- 
ically, and which includes all relevant details of kinetic 
and chemical equilibrium. The purpose of this present 
work is to better quantify whether there are observable 
consequences of the resulting reaction network. The ul- 
timate goal is to better explain the observed overproduc- 
tion (1000 times) and underproduction (3 times) of the 
abundances of 6 Li and 7 Li. 

In Section[2]we outline some of the details of the calcu- 
lations we have performed to estimate the thermonuclear 
reaction rates, and to describe the evolution and abun- 
dance of the different nuclei, X~ particles, and X-nuclei 
during the BBN. 

In Section [3J we show the results of BBN with the re- 
actions involving the nuclides that could be fo rmed with 
the em bedded X~ particles. As was done in Dedamzikl 
( 2008a|), we identify the parameter space for solving both 
the the 6 Li and 7 Li BBN abundance proble ms. The 
present work, however, differs from that of IJedamzikl 
(2008a}) in several important ways. Besides differences in 
some of the relevant reaction rates, there are two major 
difference s . One is the inclusion of X-transfer reactions. 
IJedamzikl ()2008al ) found that X-transfer reactions involv- 
ing the charge Z = 1 X-nuclei could possibly change 
the light-element abundances. That result, however, re- 
lied on reaction rates which were calculated within the 
framework of the Born approximation. Recently, how- 
ever, a deta iled study of one o f the relevant X-transfer 
reactions by IHamaguchi et al.l 1)20071 ) shows that a more 
realistic calculation gives a much smaller (factor of 10) 
X-transfer reaction rate. Based upon that study, we con- 
clude that X-transfer reactions can be neglected in this 
work. Another difference is the assumption that hadronic 
and electromagnetic (EM) decays are unimportant. That 
is, we only consider the X-nuclear reactions and assume 
that the effect the of decay products is negligible. This is 
possible if the decay lifetimes are long and/or the mass 
difference between the X" and daughter parti cle is small. 
In sup esymmetric scenarios for example (e.g. iFeng et al.l 
l2003alTbl ) . nearly equal masses would also imply a larger 
lifetime. Very large lifetimes are ruled out but lifetimes 
of or der a year are very interesting for structure forma- 
tion (|Sigurdson fc Kamionkowskll2003h . 

In Section HJ the effects X~ decay are discussed. Our 
conclusions are summarized in Section O 

2. MODEL 

In order to perform the nucleosynthesis calculations, 
we have added the relevant X-n uclei and thei r reac- 
tions to the BBN network code (jKawanol fl992l) . The 
X-nucleus Hex is particularly important for the present 
work, although the other X-nuclei are included and can 
significantly affect the results. Both proton and neutron 
captures involving the X~ particles were included when 
energetically allowed, as were "transfer reactions." We 
modified most of the thermonuclear reaction rates on the 
X-nuclei from the original rates (without X-nuclei) . The 
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two dominant effects were the lowered Coulomb barriers 
resulting from the X~ in the nucleus, and the modi- 
fied reduced mass. However, as noted below, there are a 
number of reactions for which careful additional consid- 
eration is required in order to estimate reliable reaction 
rates. 

2.1. Properties of the X~ Particle 

The X" particle is assumed to be leptonic, since there 
is motivation fr om particle physics for the existence of 
such a particle (|Cvburt et al.l 120061 ) . Identifying these 
particles as the supersymmetric partners of normal lep- 
tons implies that they would also be leptonic but of spin 
0. Such sleptonic X~ particles will be initially produced 
in pairs with X + particles. Ultimately, their annihila- 
tions in the early universe will freezeout. The final re- 
quired BBN abundance of the residual X + -X~ pairs will 
then constrain their annihilation cross section. We show 
below that this cross section is consistent with a weakly 
interacting sleptonic particle. We note, that the X + par- 
ticles, though also present during BBN will have negli- 
gible interactions with ambient nuclei (compared to X~ 
particles) due to their Coulomb repulsion and low asso- 
ciated reaction rate. Hence, they can be neglected in the 
present analysis. It is possible, however, that they could 
affect the final results to the production of electromag- 
netic and/or hadronic showers when they decay. That 
complication will be addressed in a future paper. The 
focus of the present work, however, is only the effects of 
X~ particles on BBN. 

Ultimately, the product of the X~ abundance and 
mass must be consistent with the WMAP CMB power 
spectrum. The lifetime of the X~ also has a lower 
limit because it must live long enough to allow at least 
some fraction of the X~ produced initially in the big 
bang to exist through the epoch of BBN. If the X~ 
particles decay into leptons and photons at some later 
stage, they are expected to destroy some fraction of 
the nuclei to which they had become bound during 
BBN. However, that fraction is not expected t o be 
large (IKapling hat k Raiaramanll2006t lRosenlll975l ) and 
can be neglected. 

2.2. Nuclear Binding Energies 

The reaction rates of the X-nuclei are strongly affected 
by their binding energies. In our calculations both the 
binding energies and the eigenstate wave functions of X~ 
particles were computed by taking into account the mod- 
ified Coulomb interaction with the nucleus. For this pur- 
pose, we assumed that the charge distribution of nuclides 
is Gaussian. We then solved the two-body Shrodinger 
equation by a v ariational calculatio n (Gaussian expan- 
sion method, i.e. lHivama et a l. 2003), and obtained bind- 
ing energies. The obtained values are listed in Table [TJ 
The adopted root mean square charge radii are listed 
in the second column. When experimental data exist, 
we used charge radii determined from experiment. Note 
that since the X~ particles can bound only electromag- 
netically to nuclei, their binding energies are typically 
small (~ 0.1 — 1 MeV), and are largest for heavy nuclei. 
Hence, they are not appreciably bound to nuclei until the 
temperature becomes low enough. 

2.3. Reaction Rates 



The leading term in the expression for the thermonu- 
clear reaction rates (TRR) (av) for nonresonant reac- 
tions involving nuclei with embedded X~ particles can 
be written (e.g. lBovdll2007f) as 

Mnr = (£)%bT)(0 

xE- 3/2 S(E )T 2 exp(- T ), (1) 

where E = 1.22(z^Zf /xTf) 1 / 3 keV is the energy at the 
peak of the Gamow window, S(Eq) is the " astrophysical 
5- factor" at Eq, Ub is the Boltzmann constant, T§ is the 
temperature in units of 10 6 K, and 

The astrophysical S factor contains the nuclear matrix 
element for the reaction. For some of the rates we 
take the nuclear matrix element within S(Eq) to be the 
same as for the reactions of the corresponding normal 
nuclei (NACRE; iCauehlan k Fowled Il988t iSmith et all 
Il993h . For others we can adopt a scaling relation as de- 
scribed below. Hence, the dominant corrections for the 
TRR in the above equation arises from the terms involv- 
ing the reduced mass fj, (in atomic mass units), and z\ 
and Z 2 (the proton numbers for the projectile normal 
nucleus and the target A-nucleus) . We have assumed Z 2 
to be the net charge of the bare nucleus and any embed- 
ded X~ particles, i.e. Z 2 — z 2 — n for single (n = 1) 
A~-bound nuclei, where z 2 is the atomic number of the 
bare nucleus. As noted above, we take the spin of the 
X~ particles to be zero. 

At the epoch during BBN at which the A~-nuclei 
undergo nucleosynthesis the neutron abundance is ex- 
tremely small. Nevertheless, for completeness we have 
included neutron radiative capture reactions in the re- 
action network. We adopted the known values for the 
normal nuclei whenever possible, but those values are 
not always available for the unstable nuclei of interest 
3 < z 2 < 5. For the Bex, 6^(71,7) reactions we as- 
sumed a mean value obtained f rom the 6 Li(n,7) 7 Li and 
7 Li(n ,7) 8 Li reactions (see, e.g. iHeil et al.lfl998t iBarkerl 
I1980L and references therein) of 40 /Lib at 25 keV. 

2.3.1. X~ Transfer Reactions to Produce 6 Li, 7 Li and 7 Be 

Reactions in which an X~ particle can be transferred 
can be very important in circumventing some of the reac - 
tions that would normally be inhibited (|Pospelov| [2007') . 
This is especially true of the 4 Hex(d,X~) 6 Li reaction. 
The rate for this reaction could be orders of magnitude 
larger than that of the 4 He(d, r y) e Li reaction. This lat- 
ter reaction is the main process by which 6 Li is made in 
BBN in the absence of X~ particle. Normally, however, 
this reaction is suppressed because it is dominated by an 
el ectric quadrupole t ransit ion. 

Hamagu chi et alj (|2007T ) have recently carried out a 
new theoretical calculation of the cross section for the 
X~ transfer reaction 4 Hex (d,X~) 6 Li in the context of a 
quantum three-body model. Their val ue was about a n 
order of magnitude smaller than that of lPospelovl (|2007f > . 
This difference can be traced to an exact treatment of the 
quantum tunneling in the fusion process and the use of a 
better nuclear potential. We have therefore adopted the 
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rate of iHamaguchi et ail (|2007f) as the most reliable esti- 
mate for this ra te and a ssume a factor of 3 uncertainty. 

ICvburt et al.l (|2006l ) estimated astrophysical S- 
factors for various transfer reactions including the 
4 He x (t,X-) 7 Li, 4 He x ( 3 He,X-) 7 Be, 6 U x (p,X-) 7 B e 
reactions by applying a scaling relation (|Pospelovl [2007). 



S x /S 1 ocp/a /(w 7 a ) 



2A+1 



(3) 



where Sx and 5 7 are the ^-factors for the X~ transfer 
and normal radiative processes, respectively. The quan- 
tity do is the X~ Bohr radius of 4 Hex or 6 Lix , while 
Pf is the linear momentum of the outgoing 7 Li or 7 Be. 
The quantity w 7 is the energy of the emitted photon 
of multipole order A (A = 1 for electric dipole) in the 
radiative-capture reactions. 

In the present work, however, we consider more details 
of the reaction dynamics in order to better clarify the 
differences in these reactions. First, we note that 4 He, 
6 ' 7 Li, and 7 Be occupy an s-wave orbit around the X~ 
particle (assuming the X~ particle to be much heavier 
than these nuclei). At the same time, the 6 Li nucleus is 
an a + d cluster system in a relative s-wave orbit, while 
the A = 7 nuclei are a + t and a+ 3 He cluster systems 
in relative p-wave orbits. This difference in the orbital 
angular momentum will produce a critical difference in 
the reaction dynamics between the 4 Bex(d,X~) e Li re- 
action and the 4 He x (i,^~) 7 Li, 4 Hex( 3 He,X") 7 Be, and 
G Lix(p,X~~) 7 Be reactions. In particular, the latter three 
X~ transfer reactions to produce 7 Li and 7 Be must in- 
volve a AI — 1 angular momentum transfer. This leads 
to a large hindrance of the overlap matrix element of the 
nuclear potential for the X~ transfer processes. In the 
latter three reactions the outgoing 7 Li and 7 Be in the fi- 
nal state must occupy a scattering p-wave orbit from the 
X~ particle in order to conserve total angular momen- 
tum. Thus, a realistic quantum mechanical calculation 
would de duce much smaller S x -factors than those esti- 
mated by ICvburt et al.l (|2006f ). In this article, therefore, 
we have assumed that the above three reaction processes 
are negligible. 

2.3.2. 7 Bex +p Resonant Reaction 

iBird et all (|2007h have recently suggested that the 
7 Bex(p,j) s Bx resonant reaction could occur through an 
atomic excited state of 8 Bx with a threshold energy of 
167 keV. This channel does destroy a significant amount 
of Bex and is included in our present study. 

In the previous study (Kusaka be et al.l[2007| ), we have 
suggested that a reaction channel through the 1 + , E* — 
0.770±0.010 MeV nuclear excited state of 8 B via 7 Be x +p 
-^ 8 B*(1+, 0.770 MeV) x ^ 8 B x +7 could also destroy 
Bex- However, we found that this channel is unim- 
portant fr om the estimate of bi nding energies in this 
study. In iKusakabe et ail (|2007l ). the binding energies 
and the eigenstate wave functions of the X-nuclei were 
calculated assuming uniform finite-size charge distribu- 
tio ns of radii rp = 1.2A 1 / 3 fm for nuclear mass number 
A (|Cahn & Glashowill98lh . Since the assumed charge 
radii were smaller than the measured charge radii, the es- 
timated binding energies were higher than those of this 
study. Although the nuclear resonance channel is not 
important, we show the effect of this channel on Bex 
destruction. 



The TRR for resonant radiative capture reactions is 
given approximately by 



{av) R = h 2 



2tt 
fik B T 



3/2 



L0 7 eM-E/kBT), (4) 



21 + 1 



r z r 



(2h + 1) (2J 2 + 1) r t 



(5) 



where h is Planck's constant, E is the resonance energy, 
and Ji, I2, and / are the spins of the projectile, the tar- 
get, and the resonance. The quantity Ti is the particle 
width for decay into two charged particles of relative an- 
gular momentum Z, while T 7 is the gamma decay width, 
and r to t is the total width of t he resonant state. 
Yi is approximately written (|Bovdll20"07t ) as 



3h 



r,« — 



1/2 



R \AM Uj 
x exp[6£T 1/2 



1.05 (AR Zl Z 2 



,1/2 



-7.62 (l + D (ARz^)- 1 / 2 }, 



(G) 



where R — 1A(AY 3 + A^ 3 ) x 10~ 13 cm is the interaction 
radius with Ai and Ai the atomic weights of the interact- 
ing particles. The symbol A denotes the reduced atomic 
weight A = AiA2/(Ai +A2), and M u is the atomic mass 
unit, while Of is the dimensionless reduced width. The 
height of the Coulomb barrier is 

E c = 1.44(MeV fm) 



ziZ 2 
R 



(7) 



while the usual Sommerfeld parameter b is given by 
b = 31.28z 1 Z 2 A 1 / 2 (keV 1/2 ). Correcting the charge, re- 
duced mass, and energy for the nuclides with an em- 
bedded X~ particle, we obtain the partial width of the 
8 B*(0.77 MeV)x -+ 7 Bc x +p reaction, 



T hX ~ 1.7 x 10 6 eVexp ( - 



93.9 

(Sth/keV) 1 ^ 



(8) 



where E^h is the energy level of the nuclear resonant state 
with respect to the 7 Bex plus p exit channel. 

From the conjugate analog state of 8 Li (1 + , 0.9809 
MeV) we deduce that the 8 B (0.770 MeV) first excited 
state has spin and parity T" = 1 + and that the resonant 
reaction proceeds through a p-wave (1 = 1). In Eq. © 
we also adopted 9\ = 0.82 from the TRR of this reso- 
nance used in standard BBN. Applying this proton wid th 
[Eq. ©] and T 7 = 25±4 meV (|Aizenberg-Selovell988D to 
Eqs. ([U and (J&J), we can estimate the resonant TRR for 
7 Be x +P -> 8 B*(1+, 0.770 MeV)* ^ 8 B x +7 as a func- 
tion of E. When E s» 0, T\ t x vanishes, thus oj 1 — » 
and the TRR also vanishes. As the resonance energy E 
increases, T^x also gradually increases and eventually 
becomes larger than T 7 . Thus, lo 1 converges to a con- 
stant value w 7 = (27 + l)/((2ii + l)(2/ 2 + l))r 7 inEq. ©. 
However, the exponential factor exp(—E/kBT) in Eq. (jj) 
strongly regulates the TRR for larger E. In this man- 
ner there is a most effective resonance energy E e g which 
maximizes the TRR. We found E c g w 30 keV numeri- 
cally. Hence, this resonant reaction could be an impor- 
tant, possibly dominant, new means to destroy 7 Bex- 
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The present study, however, indicates that is not the 
case. 

In the previous study (Kusakabe ct al. 2007) the cal- 
culated binding energies of the X" particle in 7 Bex and 
8 B X were respectively 1.488 MeV and 2.121 McV. If wc 
adopt these values for the energy levels of the nuclear 
excited states of 8 Bx, this 1 + state of 8 Bx becomes lo- 
cated near the particle threshold for the 7 Bex+p chan- 
nel. In this case the 7 Bex(p,j) 8 Bx reaction can pro- 
ceed through a zero-energy resonance of 8 B X at E w 
MeV, where E is the center-of-mass energy between 
the 7 Bex and the proton in the entrance channel. We 
note, however, that the estimated binding energies of 
X-nuclei, depend upon several assumptions such as the 
adopt ed charge distribution of the X ~ embedded nu- 
cleus (|Pospelovll2007HBird et alll2007h . The A=6, 7, and 
8 nuclear systems are typical clustering nuclei for which 
even a small change of the relative wave function between 
the composite nuclei can affect significantly the radiative 
capture cross sections at astrophysical en ergies as well 
as their static elec tromagnetic properties (|Kaiinolll986b 
iKaiino et al.| [T988). We therefore adopt a large uncer- 
tainty in the 1 + resona nce energy, E, from th e 7 Bex+P 
separation threshold in iKusakabe et al.l (|2007[ ) . 

To check the nuclear flow to the higher mass re- 
gion we added the 8 Bex+P -^ 9 B3f -^ 9 Bx+7 reaction, 
where 9 B™ indicates that the reaction occurs through 
an atomic excited state of 9 B. However this reaction was 
now found to be unimportant because its threshold en- 
ergy is relatively large (see Table [2]). 



2.4. Reaction Network 

We show in Fig. [1] the reaction network involving the 
nuclides with an embedded X~ particle relevant to the 
present paper. Nuclear reactions denoted by the thick 
solid arrows were found to be especially important for 
the production or destruction of the A=6 and 7 nuclides 
both i n the literature and in our present study. iPospelovl 
(|2007f ) originally proposed the main production process 
for 6 Li to be the 4 Bex (d,X~) 6 Li reaction. Similar X~ 
transfer reactions 4 Hcy (£,X~) 7 Li, 4 Hey ( 3 He,X~) 7 Be , 
6 Li x (p,X-) 7 Be were discussed by ICvburt et al.l (|2006ft 
as possible additional production processes of A=7 nu- 
clides. Here wc additionally considered the destruc- 
tion processes 6 Lix(p, 3 Hc) 4 Hcx, 7 Lix{p,a) 4 Rex , and 
7 Bex(d,pa) 4 Kex- We find, however, that these latter 
three reactions produce small effect on the BBN results 
and, as noted above, the latter three X~ transfer reac- 
tions would have little effect on BBN results when Al = 1 
h indrance is t aken into account. 

iBird et ail (|2007l ) have proposed the likely important 
destruction process 7 Bex+P -^ 8 B* X -^ 8 Bx+7, where 
8 B3f is an atomic excited state through which the ra- 
diative capture occurs. They also proposed a charged 
weak-boson exchange reaction 7 Bex -^ 7 Li+X° followed 
by 7 Li(p,a) 4 He and 7 Li(X~ ,j) 7 Lix(p,a) 4 Rex to even- 
tually destroy the A=7 nuclides. Other reactions can 
also contribute to synthesize or destroy e Li n x and 7 Li n x, 
7 Be„x (n=0 or 1) and the heavier nuclei. Our network 
code includes many reactions with nuclei up to Carbon 
isotopes. Table [2] summarizes our adopted nuclear reac- 
tion rates. 



2.5. Kinetic Rate Equations 

When dealing with the kinetic and chemical equilib- 
rium associated with X~ particles, it is necessary to 
consider the thermodynamics associated with the bind- 
ing of X~ particles. This is because it will be im- 
portant to know precisely when during BBN they be- 
come bound to nuclei, and what th eir distribution over 
the B BN nuclei is expected to be (Koh ri fc Takayamal 
2007;. We have put both recombination and ionization 
processes of X~ particles into our BBN network code 
and have dyna mically solved the a s sociat ed set of rate 
equations [as in lKohri fc Takavamal (|2007f )] to find when 
the X-nuclei decoupled from the cosmic expansion. De- 
noting a specific isotope of an element by (N,Z), its 
abundance as n(N,Z), and the corresponding quanti- 
ties for the isotopes with an embedded X~ as (N, Z)x 
and n{N,Z)x, the capture-reionization expressions for 
(N, Z) and (TV, Z)x in an expanding universe are 



dn(N, Z) 
dt 



3Hn(N, Z) 



~dn{N, Z)' 




~dn(N, Z)' 




dt 


creation 


dt 


destruction 



dn(N, Z) : 



dt 



(9) 



capture 



and 

dn(N, z) x 
dt 



3Hn(N, Z) x 



dn(N, Z) x 




'dn(N, Z) x 


dt 


creation 


dt 



+ 



dn(N, Z) x 



dt 



destruction 

(10) 



capture 



Here, the subscript "creation" refers to nuclear reactions 
that make (N, Z) and "destruction" to nuclear reactions 
that destroy it, including /3-decay for unstable nuclides. 
The quantity H is the Hubble expansion rate. De- 
tailed balance of the X-capture reac tion (7V, Z) + X~ — > 
(Nj Z )x + 7 and its inverse permits (|Kohri fc Takavamal 
l2007f ) writing the capture process as 



dn(N, Z) 



x 



dt 



capture 



(a r v) (nxn(N, Z) ~ n(N 7 Z)xn 1 ) 

. ( u ) 

where cr r is the recombination cross section, nx is the 
abundance of X~ particles, and n 7 is the number of pho- 
tons in excess of -E^ind, the X~ binding energy to (TV, Z), 
i.e., 

^ =n H^MM^H cxp r^J' (12) 

(13) 



and 



Note that since -Ebind is small, the equilibrium will favor 
unbound X~ particles until low T. Since the mass of 
the X~ particle is assumed to be > 50 GeV, the reduced 
mass for the X~ + A(N, Z) system can be approximated 
as /ix = TTiAmx / (mA + fnx) ~ tua,- This leads to 



() 



the following TRR for the first reco mbination process 
A{X-,i)Ax (jKohri fc Takavamdl27)07l ) 

2^aZ 2 (2tt) 1/2 gbind 
^ V)X ~ 3exp(4.0) M |( MX T) 1/2 

oc Z 2 E hmd m A 2 - 5 ~ Z 3 m A 2 5 . (14) 

where a is the fine structure constant. This expression is 
almost independent of mx but increases with Z. We do, 
however, obtain a different mass-dependence in the TRR 
for the second recombination process Ax(X~ ,"f)Axx ■ 

2 9 ™(Z-1) 2 (2tt) 1/2 E hind 
{<TlV)XX ~ 3exp(4.0) 

oc (Z - l) 2 S bind m^ 2 - 5 ~ (Z - l) 3 m x 2 - 5 . 

(15) 

Here we obtain a mass dependence because fixx = 
nT-Axm-x /(mAX + mx) ~ mj/2. Since is typically 
much larger than the mass of the light nuclei mx 3> tua, 
the rate for the second or higher-order recombination 
process is hindered. 

3. RESULTS 
3.1. BBN Calculation Result 

Figures and^) illustrate the results of a BBN cal- 
culation in which the X~ abundance is taken to be 10% 
of the total baryon number, i.e. Yx = nx/n^ = 0.1, 
where nx is the number density of the X~ particles and 
rib is the averaged universal baryon-number density. Re- 
sults shown in Fig. [2^ are for normal nuclei while Fig. [2jz> 
sh ows those for the X-nucl ei, Note that we adopt the rate 
of lHamaguchi et all (|2007l ) for 4 HexKX"), and the re- 
actions 4 Be x {t,X-), 4 He x ( 3 He,X-), and 6 U x (p,X-) 
are taken to be negligible. 

The abundances for the normal nuclei 6 Li, 7 Li and 7 Bc 
in the interval Tg ~ 0.5 — 0.2 are seen to be close to their 
usual BBN values. This is because at higher temper- 
atures the nuclear statistical equilibrium does not par- 
ticularly favor the production of the weakly bound X- 
nuclci. Thus, their effect on the abundance is minimal. 
There are some changes, however, once the X-nuclei ap- 
pear. Because of the larger binding energies, the X~ 
particles bind first to the heaviest nuclides, such as 7 Li 
and 7 Be produced in normal BBN (first recombination 
of the X~ particles). These recombinations occur at 
around T g = 0.3 (for 7 Li) and T 9 = 0.5 (for 7 Be), re- 
spectively. An increase in the Bex abundance by the 
recombination can clearly be seen in Fig. [2Jd. Some- 
what later, at around Tg — 0.1, the X~ particles are 
captured onto 4 He, as can be seen in Fig. Then 
a new round of nucleosynthesis of the X-nuclei, involv- 
ing the reaction 4 Hex(rf,X~), produces normal nuclei 
6 Li (in Fig. Hi) as well as 6 Li x (in Fig. gb). How- 
ever, the most notable feature of these results is that 
the 6 Li is not easily destroyed by the 6 Li(p,a) 3 He reac- 
tion, which destroys nearly all of the 6 Li produced in 
standard BBN. This is because the X~ transfer reac- 
tion restores the charge of the normal nucleus 6 Li, which 
has a Coulomb barrier which is too high at these tem- 
peratures for its destruction. Thus, the large abundanc e 
ratio of mass 6 to mass 7 is preserved ()Pospelovl 120071 ) . 



In addition, however, Be is destroyed by the nuclear re- 
actions that occur after the recombination (Tg ~ 0.3), 
mainly 7 Be(X~,7) 7 Bex(p,7) 8 Bx- This is explained in 
detail below. 

The calculated BBN abundances of the mass 6 and 7 
nuclides and other light nuclei depend strongly on the 
X~ abundance. The above discussion applies only to 
the case of relatively abundant X~ particles. In order 
to study the sensitivity of the 6 Li and 7 Li (= 7 Li+ 7 Be) 
abundances to the X~ abundance, nx, we carried out a 
series of BBN calculations in which nx was varied over 
a wide range. 

Figure [3] shows the calculated abundances of the mass 
6 and 7 nuclides. These are plotted as 6 Li/H (solid curves 
of positive slope) and 7 Li/H (horizontal solid curve), as a 
function of the initial X~ number fraction fx relative to 
the cold dark matter (CDM) abundance, which we define 
as 

, Y x r^/50GeV\ 
0.09 n C DM V "»cdm / 

This definition follows if we take number density of CDM 
particles as that inferred from the CDM closure content 
(fi CDM = 0.2) deduced from the WMAP data. For this 
closure parameter, the CDM number density, for any 
value of mass jticdm of the CDM particle, can be written 
ncDM = 0.09 rib X (50 GeV/mcDivi)- So, our definition 
of fx is for a fiducial CDM mass of 50 GeV and is easily 
scalable to other CDM masses. 

In the calculations shown in Fig. [3] we have assumed 
that the X" particle has a mean lifetime much longer 
than the typical time scale for BBN in the presence of 
X~ particles, i.e tx 3> 5 hours. Below we also con- 
sider the case of a mean lifetime which is shorter than 
5 hours. In this figure we have taken into account the 
theoretical uncertainty in the X~ transfer reaction cross 
sections. For 6 Li the upper curves correspond to the 
yields with our assumed rate fo r the 4 Hex(d 7 X _ ) 6 Li re- 
action (|Hamaguchi et aill2007f ) multiplied by a factor of 
3, and the lower curves correspond to that same rate 
divided by a factor of 3. The middle curve for 6 Li corre- 
sp onds to our best guess for that crucial rate, i. e., the rate 
of lHamaguchi et all (|2007l ). As noted above. iBird et al.l 
(2007) suggested that the recombination of 7 Be and X" 
together with the 7 Bex (p, 7) 8 Bx reaction taking account 
of the resonant contribution from both reactions would 
destroy a considerable amount of 7 Be; this process was 
included in all of our calculations. 

It is clear from Fig. [3] that the 6 Li abundance in- 
creases monotonically with increasing fx- This is a 
consequence of the fact that 6 Li is mainly pro duced by 
the 4 H ey(ri,X~) 6 Li reaction as proposed in iPospelovl 
(2007) for almost entire range of /x-values. For 
fx 10~ 9 , however, the standard BBN processes 
4 He(rf,7) 6 Li and 3 He(t,7) 6 Li are the main reactions t o 
make 6 Li (jFukugita fc Kaiino|[l990l : ISmith et al J 1 19931) . 
In the region of fx ^ 0.1, however, a departure from 
the linear increase due to 6 Li(X _ ,7) 6 Lix(p, 3 He) 4 Hex is 
observed. 

During normal BBN, at the baryon to photon ratio 
deduced from CMB measurement by WMAP, 7 Li is pro- 
duced mainly as 7 Be. We find that the 7 Be (and 7 Li) pro- 
duced during the standard BBN epoch captures X" par- 
ticles to form 7 Bex (and 7 Lix) during the first recombi- 
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nation. The 7 Be is then destroyed by the 7 Bex(p,'~f) s Bx 
reaction. In the parameter region of lower X~ abundance 
fx < 0.1, 7 Be and 7 Li are produced in the standard 
BBN by the 4 He(i,7) 7 Li and 4 He( 3 He,7) 7 Be reactions. 
They are not, however, destroyed and remain almost un- 
changed. 

3.2. Constraints on the Primordial 6 Li Abundance 

Constraints on the primordial 6 Li abundance must 
be inferred from the observed plateau as a function of 
metallicity on MPHSs. The actual primordial abundance 
of 6 Li could be hig her than the recently detected high 
plateau-abundance (|Asplund et al.l [2006) of 6 Li (lower 
horizontal dashed curve in Fig. [3]). This is because 
stellar processing could have deplete d an initial surface 
abun dance. It is expected in models (|Pinsonneault et"aLl 
120021 ) of stellar structure and evolution that both 6 Li and 
7 Li abundances decrease because materials on the stellar 
surface might be convected to regions of sufficiently high 
temperature that the fragile 6 Li and 7 Li are par tially de- 
stroyed (e.g. lLamberll2004lRlc~hard et alJl2005f ). There- 
fore, the plateau level for the observed abundances of 
6 Li/H and 7 Li/H in MPHSs should be considered a lower 
limit to the primordial abundance. For this lower limit 
we take the 3 a lower limit to the mean plateau value 
times a factor of 1/3 for systematic uncertainties giving 
6 Li/H > 1.7 x 1CT 12 . We include this additional factor 
because there may be additional systematic uncertain- 
ties due to the sensitivity of the inf erred 6 Li abundanc e 
to the model atmosphere employed ()Cavrel et al.ll2007t ). 

The upper limit is more difficult to estimate. Be- 
cause 6 Li could be more easily destroyed in stars than 
7 Li, its upper limit should be higher than the upper 
limit to the 7 Li abundance. Even so, there are con- 
straints on the degree of stellar processing for both 7 Li 
and 6 Li from the limits on the dispersion of the plateau. 
A large degree of stellar destruction would be sensitive 
to the varying degrees of meri dional circulation in the 
stars (|Pinsonneault et all l2002t ) and hence should pro- 
duce a large dispersion in observed abundances. How- 
ever, the fact that the observed dispersion in 7 Li is 
greater than that observed for 6 Li (albeit on a limited 
data set) s uggests that the destruct ion may not be very 
significant (|Pinsonneault et al. 2002) . For this reason we 
adopt a conservative upper limit of a factor of 10 above 
the mean plateau value giving 6 Li/H < 7.1 x 10 -11 . 

3.3. Observational Constraints on the X~ Abundance 
3.3.1. Case of Longer Mean Life of X~ Decay 

We consider several constraints on the X~ abundance 
from the observed Li isotopic abundances. We first dis- 
cuss the case of a longer mean life for the X~ particle 
(tx > 5 hours). In this case, the yields of BBN are not 
strongly affected by the value of Tx- Our adopted lower 
limit to the 6 Li/H ratio is then satisfied by the calculated 
abundance of 6 Li/H (thick solid curve in Figs. [3]and[2]) 
for / x >2x 10- 6 . 

In order to include the uncertain depletion effect of 
fragile lithium in stellar atmospheres, we display in Fig. 
131 the calculated primordial 6 Li and 7 Li abundances rela- 
tive to the mean plateau levels in MPHSs, i.e. (^Li/H)/ 
(' 4 Li/H)MPHS- The gray regions enclose our adopted un- 
certainty due to the X~ transfer reaction cross sections 



as in Fig. [3] Here and below, we take, as our recom- 
men ded network, the reactio n rate of 4 Hex(d,X~) 6 Li 
from lHamaguchi et all (|2007ft . 

The parameter regions of initial X~ abundance, 2 x 
10~ 6 </x<lx 10~ 4 for 6 Li and fx < 1 for 7 Li, respec- 
tively, are plausible regions in this context. In addition, 
since 6 Li is more easily destroyed in proton burning than 
7 Li, the primordial abundances should satisfy the in- 

equality ( 7 Li/H)/( 7 Li/H) M p H s < ( 6 Li/H)/( 6 Li/H) M p H s- 
From Fig. 3] we thus find the concordant parameter re- 
gion of 

4 x 10~ 5 </ x <9x 10~ 5 , 
(i.e. 4 x 10~ 6 < Y x < 8 X 10~ 6 ), (17) 

which is bounded by the vertical solid line on Fig. |4j 

Here, we used Eq. (| to convert the inferred lim- 
its on fx into the limits for Yx- We then deduce from 
Eq. (fT7|) and Fig.|4]a possible depletion factor, d( A Li), of 
the primordial abundances of 6 Li and 7 Li in Population 
II metal-poor halo stars of, 

d( 6 Li) < 10, (18) 
d( 7 Li) < 4, (19) 

where the upper limit to 6 Li depletion comes from our 
adopted upper limit to the primordial 6 Li abundance. 
We note that these depletion factors are consistent with 
the known thermonuclear reaction rates for p+ 6 hi and 
p+ 7 Li, assuming partial depiction of lithium in the so- 
lar atmosphere; the former is about a factor of 80 larger 
than the latter. We also point out that in our present 
scenario for BBN including the X~ particles, their abun- 
dance parameter region, Eq. (fT7j) can be a solution to the 
Li abundance discrepancy between the standard BBN 
prediction and the observations of MPHSs, while still 
satisfying the independent abundance constraint on the 
primordial 7 Li abundance. It is to be noted, however, 
that there still remains a possible controversy that the 
depletion factor of G?( 7 Li) = 3 — 4 may be too large to ac- 
commodate the observed small dispe rsion in the plateau 
abun dance level detected in MPHSs (iPinsonneault et all 
l2002h . 

3.3.2. Case of Shorter Mean Life of X~ Decay 

To study the effects of X~ decay we introduce the life- 
time Tx of the X~ particle. Then there are two param- 
eters Yx and Tx ■ Even in the case when the mean life of 
the X~ particle is nearly equal to or slightly shorter than 
five hours (tx < 2x 10 s) the recombination of 7 Be and 
X~ particles still enriches the 7 Bex abundance (Fig. [2ja) . 
This is because the recombination of 7 Be and X~ parti- 
cles occurs at earlier times when the cosmic temperature 
is Tg ~ 0.5. However, in this case the recombination 
of 4 He cannot produce abundant 4 Hex- This is because 
the recombination of He occurs at at a lower tempera- 
ture Tg ~ 0.1 when the cosmic time is of the same order 
of tx ~ 10 4 s and is sensitive to the decay. Thus, the 
6 Li production is reduced because the 4 Hex(d,X~) 6 Li 
reaction is strongly hindered. 

In this case, however, th e resonant reacti on process 
7 Be*+p ^ 8 B*jl ^ 8 B x +7 (|Bird et al.|[20TJ7T) can more 
effectively destroy 7 Bex- We assumed that 8 Bj inter- 
converts to s Bex by /3-decay with a rate of 8 B /3-decay 



multiplied by the correction term (Qx /Q) 5 j where Q and 
Qx are the Q-values of standard /3-de cay and that of 0- 
decay for X-nuclei. We adopt here the lHamaguchi et al.1 
(2007) rate for 4 He x {d,X-) 6 U. 

In Fig. [H the contours of c?( 6 Li) [solid curves] and 
c£( 7 Li) [dashed curves] are shown. The upper and lower 
solid curves correspond to the abundance level which 
satisfies our adopted constraint (1.7 x 10~ 12 < 6 Li/H 
< 7.1 x 10 -11 ) as discussed in subsec. 13.21 from the 
abund ance of MPHSs, ( 6 L i/H) M PHs = (7.1 ± 0.7) x 
10~ 12 ([Asplund et al.l 12006). Thus the upper right re- 
gion of the figure is excluded for 6 Li overproduction. The 
right-upper side from the lower solid curve indicates the 
region for which the 6 Li abundance is higher than that 
observed in MPHSs. 

The three dashed curves correspond to d( 7 Li)=1.55, 2, 
3, from right to left, respectively. In Fig. [5] there is no 
7 Li over-destruction region where 7 Li abundance is below 
the observed mean value. <i( 7 Li)=1.55 cor responds to 
the 1 sigma upper limit for MPHSs value (jRvan et al.l 
2000). In the right -upper side from the dashed curve 
for d( 7 Li)=2, the resulting 7 Li abundance is lower than 
7 Li/H« 2.5 x 10~ 10 . Therefore, we conclude that BBN 
with negatively charged particles provides a simultaneous 
solution to the 7 Li overproduction problem and the 6 Li 
under production problem, [as also deduced bv lBird et al.l 
(|2007h ] in the parameter region 

Y x > 0.9, t x ~ (1.0- 1.8) x 10 3 s. (20) 

For tx > 10 5 s and Yx > 3, the calculated abun- 
dance of 7 Li increases slightly. In this region the 
6 Lix (p, 7) 7 Bex reaction produces some amount of 7 Bex . 
However, this parameter region is uninteresting due to an 
extreme overproduction of 6 Li compared with observed 
abundances. 

In Fig. [5] there is a solid line corresponding to 
d( 6 Li)=4. The condition 

( 7 Li/H)/( 7 Li/H) MP HS < ( 6 Li/H)/( 6 Li/H) M p HS or equiv- 
alently c?( 7 Li)< <i( 6 Li) is satisfied in the gray region. The 
gray colored region in Fig. [5] [cf. Eq. (fl8|) ] is the most in- 
teresting and relevant parameter region in order to solve 
both the 6 Li and 7 Li problems. The wider vertical band 
of 4 x 10~ 6 < Y x < 8 X 10~ 6 and t x > 10 4 s corresponds 
to the solid box bounded by the vertical solid line in 
Fig. 0] and Eq. (|17p . This region of the parameter space 
solves only the 6 Li problem, but leaves the 7 Li problem 
unresolved. 

Finally, we consider another case where 7 Bex converts 
to 7 Li by a weak charged current transition from X~ 
to X°, i.e. 7 Bex -^> 7 Li+X°. This de cay quickly trans - 
forms 7 Be x to 7 Li (type II model in iBird et all 120071 ). 
In this model, the rate of recombination effectively de- 
termines the 7 Be^ 7 Li conversion rate induced by X~. 
The results of this case are shown in Fig. [5J The general 
features of Fig. [5] are very similar to Fig. [5] However, 
due to a slightly stronger destruction rate due to the 
7 Bex -^ 7 Li+V° decay followed by the 7 Li(p,a) 4 He and 
7 Li(X _ ,7) 7 Lix(p,a) 4 Hex reactions, the contours of the 
7 Li abundance are systematically shifted toward smaller 
values of Yx- The parameter region which solves both 
the 6 Li and 7 Li problems also slightly shifts to 

Y x w 0.04 -0.2, T X w (1.4 - 2.6) x 10 3 s. (21) 

4. DISCUSSION 



4.1. Constraint on Dark Matter Particles 

If the dark matter particles, which we denote Y° in 
this article, are the decay products of X~ particles, the 
cosmological parameter JIcdm inferred from the WMAP- 
CMB data can be used to constrain the mass of Y° 
when combined with our abundance constraints on the 
X~ particles from Eqs. (|20p and (|2"Tj) . We suppose 
that X~ decays to a dark- matter Y° particle and any 
residues, and Y Y = Y x . The WMAP-CMB constraint 
on f^cDM = 0.2 corresponds to Yymy < 4.5 GeV, i.e. 



Y 



x ^ 



Yy < 



4.5 GeV 

my 



(22) 



The calculated abundance constraints on Y x (Eqs. [201 
and [22] can therefore be used to constrain the mass my. 
Note, that the calculated result in the present study does 
not particularly depend on the assumed mass of the X~ . 
In fact, only the second recombination rates and the nu- 
clear reaction rates between V-nuclei depend on the mass 
of the X~ (see discussion in subsec. 12. 5p . However, the 
main production and destruction processes of 6 Li, 7 Li, 
and 7 Be are completely free from these processes. We 
can thus consider the general constraint on the number 
fraction of X~ particles, Yx- 

The most interesting solution to both the 6 Li under- 
production and Li overproduction involves the parame- 
ter space defined in Eqs. ([20]) and ([21]). Using [Eq. ([20]) ] 
in which one includes the destruction react ion process 
7 Be x +P -> 8 B™ ^ 8 B x +7 (|Bird et al.l 12001 . the mass 
of the dark matter particle Y° would be constrained to 
be 

my < 5 GeV . (23) 

However, using [Eq. ([2T])] in which one includes the 
7 Be x ^ 7 Li+V° process ([Bird et al.ll2007h the allowed 
mass range increases to 



my < 20 - 110 GeV, 



(24) 



The long lifetime for X~ decay is consistent with the 
mass of the decaying particle X~ being close to the mass 
of the daughter particle Y°. In this case one can deduce 
a constraint on the X~ mass of 



m x < O(100 GeV). 
4.2. Initial X~ Abundance at BBN Epoch 



(25) 



Here, we consider a simple estimation of the initial X~ 
abundance after the BBN epoch taking X + and X~ pair 
annihilation into account. When X + and X~ particles 
are abundant at high temperature, both pair creation 
and annihilation equilibrate. However, as the universe 
expands and cools, the annihilation process proceeds un- 
til these particles freezeout at some relic X~ (X + ) abun- 
dance. A calculation of this is analogous to the well 
known we akly-interacting massiv e particle (WIMP) cal- 
culation of lKolb fc Turner! (|1990f ). 

Following their derivation, the time evolution the X- 
to-photon ratio, r\x = nx/n~ ( = r/Yx, can be written, 



di] 



— - -n 7 (av)r] x . 



(26) 



However, unlike the iKolb fc Turner! (jl990f ) calculation, 
the X~ annihilation cross section is given by the electro- 
maganetic formation of a positronium-like bound state 
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(X~,X + ) system, not by the weak interaction associ- 
ated WIMPs. Therefore, the annihilation cross section 
does not scale the same as the WIMP annihilation rate, 
{av)wiMP {^F rn x)-i wnere Gf is the Fermi coupling 
const ant. Rather, the much larger annihilation cross sec- 
tion ([Bird et al.ll2007[ ) of interest here is given by, 



4.4. 



(av) 



2l«> 7r 3/2 a 3 



(27) 



3exp(4)m^ /2 T 1 /2 
The solution of Eq. (f2"6"|) is then 

W (T) = ^ + 2flr 1 M->.(i-(|) 1/a )) • 

(28) 

where the subscript i refers to values at some initial tem- 
perature Ij. 

Assuming that r/xi ~ 0(1) and that the baryo n-to- 
photon ratio is r\ = 6.0 x 10~ 10 (jSpergel et al.ll2007ft . the 
freeze-out abundance of X~ can be written; 



1 



1 



x r — 



r) 2H i 1 n li (av) 

0.0057 ( -JtL_) 
V10.75/ 



V 



m x 



■50 Gev) 



6.0 x 10- 10 , 

3/2 ( J! \ ~V 2 



50 



GeV/ 



(29) 



where is the total number of degrees of freedom of 
the relativistic particles. This approximation indicates 
that the charged massive particles that existed in the 
early universe might remain as relics in abundance of 
order Yx ~ 0.01. This is not very different from the Y x - 
value of the most interesting solution to the 6 Li and 7 Li 
problems in Eq. (j20|) or Eq. (j2l"j) which we found, i.e. the 
solution which leads to the simultaneous destruction of 
7 Be and production of 6 Li at levels that would produce 
concordance with observations. 

4.3. Direct Destruction of Nuclei by X~ Decay 

Having shown how primordial X~ particles could re- 
sult in modified Li isotope production, we need to con- 
sider what would happen to those nuclei which retained 
an embedded X~ particle wh en the X~ particle decays. 
This decay was assumed by iKa plinghat & Rajaraman 
(2006) to interact sufficiently strongly with the host nu- 
cleus that it would induce nucleon emission. For exam- 
ple, the decay could knock out a proton or neutron from 
4 Hex, producing either 3 He or 3 H. These nuclides could 
then interact with other 4 He nuclei to produce 6 Li. 

However, this interaction was found by those authors 
to occur too in frequently, w hich is consistent with ear- 
lier studies (e.g. lRosen|[l975l) . It was therefore concluded 
that the resulting decay products would have little ef- 
fect if they were comprised only of leptons and photons. 
Even so, they might destroy some fraction of the nuclei 
in which the dec ays occurred if one of the decay prod - 
ucts was a pion (|Rosenlll975l : iKohri fc Takavamai r2007). 
However, for the X~ decays, the decay products might 
be expected to be electroweakly interacting particles, and 
the X~ decays might not be expected to significantly af- 
fect the host nuclei. We therefore neglect decay-induced 
destruction in the present study. 



6 Li Production by Electromagnetic Energy 
Injection from X~ Decay 



If the X~ particles decay to charged leptons and any 
other residues, the high-energy leptons thus produced 
could interact with background photons to lose their en- 
ergy and induce additional nonthermal nucleosynthesis. 
In this section we argue that such nucleosynthesis (if it 
does occur) does not significantly alter the parameter 
constraints deduced above. 

There is also a viable possibility that the decaying 
product is an electron. Such a high energy electron 
quickly interacts with background photons (e^ + 7 — > 
e ± + 7 i.e. inverse Compton scattering). The newly 
produced energetic photons will also interact with back- 
ground photons (7 + 7 — * e~ + e + i.e. pair pro- 
duction). In this way high energy electrons can lose 
energy by making an electromagnetic cascade shower. 
Suc h a shower could trigger n on -thermal nucleosynthe- 
sis (iDimopoulos et al.l Il988allbl. Il989t Uedamzikl 1 2000 1 ; 
Cvburt et al.ll2003l:lJedamzikll2004l:lKawasaki et al.H2005l : 
Kusakabe et al. 200$). 

In iKusakabe et al.l (|2006[) the conditions on the pa- 
rameters of X~ particles which lead to a resolution 
of the 6 Li underproduction problem were delineated. 
These are Y x m x > Y X E BM ~ 10~ 4 - 10~ 3 GeV and 
rx ~ 10 8 — 10 12 s, where -Eem is the generated energy 
in the electromagnetic decay process. The constraint on 
Yx for a given X~ mass mx is then 



(lO" 4 - 10" 3 GeV) 
m x 



(30) 



In the present paper we have considered several de- 
struction processes of 7 Bcx- Here we note that both 
Fig.[5]and Fig. [5] indicate the same constraint Yx < 10 -5 
from X~ BBN for tx ~ 10 s — 10 12 s which is the inter- 
esting lifetime range in the previous work based upon 
electromagnetic X~ decay (Kusaka be et al.1 12006). The 
electromagnetic decay of X~ particles triggers late time 
6 Li production in amounts of 1 < <i( 6 Li)< 10 for these 
parameter regions. Therefore, the compatible condition 
which avoids an overproduction at the earlier epoch and 
still allows for non-thermal production of 6 Li at the later 
time is only d( 6 Li)< 10. From this we deduce, 



(10- 4 - lO" 3 GeV) < 



x 10' 



m x 



(31) 



which places a lower limit on the mass of the X particle 
of 



m x > 10 - 100 GeV. 



(32) 



When mx satisfies this lower limit, the electromagnetic 
decay of X~ particles can also operate at times long 
after the BBN epoch. Clearly, electromagnetic and/or 
hadronic decays may also contribute to the final com- 
puted abundances. And, in a subsequent work we pro- 
pose to investigate this. However, for now our goal has 
been to clarify the role of X~ reactions and confirm that 
this effect alone can explain the observed lithium abun- 
dances. Note that in the case where the stable daughter 
dark matter particle has nearly the same mass as the 
X~ particle, there is not enough Q-value to produce a 
hadronic or electromagnetic shower. 
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5. CONCLUSIONS 

We have investigated light-element nucleosynthesis 
during the big bang in the presence of massive, 
negatively-charged X~ particles. Such particles would 
bind to light nuclei in the early universe. As suggested by 
many authors, they would facilitate BBN by enhancing 
the nuclear reaction rates both by reducing the charge of 
the bound A-nuclei, and by enabling transfer reactions 
involving the X~ particles. We considered the recombi- 
nation processes of X~ particles and normal nuclei. Our 
conclusions are as follows: 

First, as suggested in previous studies, the X~ parti- 
cles greatly enhance the production of 6 Li. The main pro- 
duction process of 6 Li is the sequence of Hex production 
through X~ capture on 4 He followed by the X~ transfer 
reaction 4 Hex(d,X~) 6 Li. The resultant 7 Li abundance, 
however is almost the same as normal BBN value unless 
there is a large X~ abundance Yx ^ 0.1. 

Secondly, when the lifetime of the X~~ particle is 
much longer than the period of normal BBN, the 6 Li 
abundance monotonically increases with the X~ par- 
ticle abundance, except for very small X~ abundance 
(Yx < 10~ 10 ). On the other hand, the 7 Li abundance is 
nearly independent of the X~ particle abundance unless 
the X~ particle abundance is larger than ~ 0.1 times the 
total abundance of baryons. In this case the 7 Li abun- 
dance decreases with the X~ particle abundance due to 
the resonance reaction of 7 Bex(p,7) 8 Bx which reduces 
the 7 Li abundance. 

Thirdly, the 6 Li/H and 7 Li/H observed in MPHSs can 
constrain the lifetime and abundance of an X~ particle. 
These observational constraints require the lifetime and 
abundance to be in the ranges of Tx ~ (1.0 — 1.8) x 10 3 s 
and Y x > 0.9. When the reaction 7 Be x -* 7 U+X° is 
also taken into account, these ranges change to tx ~ 
(1.4 - 2.6) x 10 3 s and Y x « 0.04 - 0.2. Therefore, intro- 
ducing X~ particles with an adequate lifetime and abun- 
dance can be a solution for both of the factor of ~ 1000 
underproduction of 6 Li and the factor of 3-4 overproduc- 
tion of 7 Li in standard BBN. 

Fourthly, a constraint on the X~ particle mass can 
be made from the dark-matter content deduced from 
the WMAP analysis of the CMB. If the abundance of 
X~ particles is > 0.1 — 1 as summarized above, 
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the mass of dark matter particles Y , which are pro- 
duced from the decay of X~ particles, turns out to be 
m Y ^; 10 — 100 GeV, thus leading to the constraint 
mx < O(100 GeV) when mx ~ my . 

In summary, although several possible solutions have 
been proposed to solve the underproduction problem 
of 6 Li, they do not necessarily resolve the overproduc - 
tion problem of 7 Li sim ultaneously. iPospelovl (|2007f ) 
and ICvburt et al.l (|2006f l proposed the X~ transfer re- 
actions to pro duce 6 Li and 7 Li. Rega rding 6 Li produc- 
tion, however. Hamagu chi et all (|2007f ) have shown that 
the assumed reaction cross s ection is no t as la rge as the 
value deduced by Pospelov. iBird et all (|2007l ) proposed 
new destruction processes of 7 Bex through atomic exci- 
tations, and we here investigated yet another destruction 
process through a 8 R X nuclear excited resonance in order 
to resolve the 7 Li overproduction problem. The destruc- 
tion efficiency of these newly proposed processes depends 
on the excitation energies above the p+ 7 Bex separation 
threshold. We then found that a better estimation of 
the binding energies leads to negligible effect from the 
nuclear excited resonance. It would be useful, however, 
in future work to predict more precisely the binding en- 
ergies and excited states of exotic A-nuclei and their re- 
action cross sections utilizing a more realistic quantum 
mechanical treatment. 
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Fig. 3. — Calculated abundances of B Li/H and 7 Li/H as a function of the initial X~ abundance parameter fx defined in the text. The 
gray band bounded by curves of 6 Li show the uncertainty from the rate for X~ transfer reaction. The dashed lines indicate the mean 
values observed in MPHSs. The dashed and solid boxes indicate the range of fx consistent with our adopted limits on the abundance of 
6 Li/H observed in MPHSs. 
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Fig. 5. — Contours of constant lithium abundance relative to the value observed in MPHSs, i.e., d( 6 Li) = SLiCaUyejjObs ^ so jj(j 
curves) and d( 7 Li) = 7Li Cal 7 7 Li° bs (dashed curves). The adopted abundances are 7 Li/H= (1.23;to'!|) X lO" 10 l lRvan et al.|[2000T) and 
6 Li/H= (7.1 ± 0.7) X 1(T 12 l|Asplund et alj|2006h . Thin solid lines around the lines of d( 6 Li) = 1 curve enclose the 1 a uncertainty in the 
adopted observational constraint based upon the dispersion of the observed plateau. In the gray region, the condition d( 6 Li)> d( 7 Li) is 
satisfied. 
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Fig. 6. — Same as in Fig. [5] when the charged-current decay of 7 Bex — > 7 Li-\-X° is included. Thin dashed lines around the lines of d( 7 Li) 
= 1 curve enclose the 1 a uncertainty in the adopted observational constraint based upon the dispersion of the observed plateau. 
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TABLE 1 

Binding Energies of X" Particles to Nuclei 



nuclide 


r RMS ( fm )a 


References 


S Bind (MeV) 


lu^ 
n x 


n -t- n nrt7 




u.uzo 


IX 


9 1 1 K 4- n nnfi 

Z.llO IE U.UUO 


9 




a X 


1 . f 00 ZLZ U.UoD 


o 


n n79 

U.U I z 


3LTo„ 


i.yoy in u.uou 


*J 
o 


n 9fis 
u.zuo 


4 He x 


1.80 ± 0.04 


4 


0.343 




9 ah -\- n n 1 ? 

Z.^o ZC U.UO 


/I 
1 


U.oUD 


7 Lix 


2.43 ± 0.02 


4 


0.882 




2.42 ± 0.02 


4 


0.945 


6 Bex 


2.52 ± 0.02 b 


1 


1.234 


7 Be x 


2.52 ± 0.02 


4 


1.324 


8 Be x 


2.52 ± 0.02 b 


4 


1.401 


9 Be x 


2.50 ± 0.01 


4 


1.477 


7 Bx 


2.68 ± 0.12 c 


5 


1.752 


8 Bx 


2.68 ± 0.12 


5 


1.840 


9 B X 


2.68 ± 0.12 c 


5 


1.917 



No te. — Re fere nces: 1= lYao fc et all pOOHl : 
2= ISimon et all lfl98ll); 3=TUNL Nuclear 
Data, http: / /www . tunl .duke.e du/NuclData 4= 
ITanihata etaLl U988I ); 5= lFukuda et al.l fL999l) . 
a Root mean square charge radius 
b Taken from 7 Be radius 
c Taken from 8 B radius 
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TABLE 2 

Nuclear Reaction Rates of X-Nuclides 



Reaction 


Reaction Rate (cm 3 s 1 mole 1 ) 


Reverse Coefficient^ 


Q (MeV) 




7 9 


o. yo 


on 


Qtt / j \/1tt 

3 He x (d,p) 4 He x 


3.9 x 10 10 T 9 ~ 2 / 3 exp(-5.36/Tg^ 3 ) 


8.49 


18.428 


6 Lix(n, 7 ) 7 Li x 


5.2 x 10 3 


1.48 


7.327 


6 Li x (n,t) 4 He x 


1.7 x 10 8 


0.577 


4.321 


6 Li x (p, 7 ) 7 Be x 


5.5 x 10 5 T~ 2/3 ex P (-6.74/T 9 1/3 ) 


1.48 


6.124 


6 Li x (p, 3 He) 4 He x 


2.7 X 10 10 T 9 _2/3 exp(-6.74/T 9 1/3 ) 


0.577 


3.557 


7t ■ / \8t - 


5.2 X 10 


1.58 


2.096 


7t ■ / \8t> 


1.3 X 10 T g ' exp(— 6.74/Tg' ) 


7.89 


17.774 


7 hi x (p,a) 4 Hex 


9.0 x 10 8 Tg _2/ ' 3 exp(-6.74/T g 1 / 3 ) 


1.00 


16.808 


8 Li x (p, 7 ) 9 Be x 


1.0 x 10 5 T g J/ exp(-4.25/T 9 1/J ) 


2.47 


17.420 


6 Bex(",7) 7 Bex 


5.2 x 10 3 


0.493 


10.766 


6 Be x (n,p) 6 Li x 


2.7 x 10 9 


0.333 


4.642 




5.2 x 10 3 


7.89 


18.977 


7 Be tc (n.ui) Li v 


2.7 x 10 9 


1.00 


1.202 




I fi y 1 8 T _ 2//3 OYnf 8 Sfi /T 1 / 3 \ 
i.o a iu i g exp^ — O.OU/ ig ^ 

II ^ v in^T" 1 3/2 _ i no /TU \ 
-f-l.D X 1U lg exp^ — L.zfZ/lQ) 


± .00 




7 Be x (rf,p) 8 Be x 


8.9 x 10 11 T g _2//3 exp(-11.13/Tg /3 ) 


16.97 


16.752 


'Bex(d,P») Hex 


8.9 x 10 11 T g -2/ ' 3 exp(— 11.13/Tg' 3 ) 


2.15 


15.786 


8d„ _,\9 


O.Z X 1U 




1 7/1 1 
1.(41 


8 Bex(p,7) 9 Bx 


1.4 x 10 6 Tg" 3/2 exp(-5.60/T 9 ) 


2.47 


0.331 


Qr) / \1|)n 

J Bex(^,7) Bex 


5.2 x 10 J 


7.89 


6.913 


a Be x (p,7) B x 


1.2 x lO'Tg ' exp(— 8.83/Tg' ) 


1.13 


7.145 


1 t-» / \11t-» 


1.2 x 10 7 T g ~ 2/ ' 3 exp(-8.83/T 9 1 ^ 3 ) 


0.493 


11.750 




5.2 x 10 3 


1.13 


8.556 


D X\PH) ^X 


1 Q y 1fl 5 T~ 2 ' 3 pvnf-1fl Tfl/T 1 ' 3 ! 


7.89 


4.727 


10 B x (n,7) u B x 


5.2 x 10 3 


3.45 


11.517 


10 B x (p,7) 11 C x 


1.9 x 10 5 T~ 2/3 exp(-10.70/T Q 1/3 ) 

y r v / y / 


3.45 


9.366 


11 B x (p, 7 ) 12 Cx 


1.9 x 10 7 T Q " 2/3 exp(-10.70/T Q 1/3 ) 

y r \ / y / 


7.89 


16.638 


10 C x (n,7)"C x 


5.2 x 10 3 


0.493 


13.195 


n CxM 12 Cx 


5.2 x 10 3 


7.89 


18.789 


4 Hex(rf, 7 ) 6 Lix 


2.2 X 10 1 Tg" 2/3 exp(-5.35/T g 1/3 ) 


2.79 


1.936 


3 Hex(a, 7 ) 7 Bex 


2.9 x 10 6 Tg" 2/3 exp(-10.70/T g 1/3 ) 


3.95 


2.642 


4 Hex(i, 7 ) 7 Lix 


2.0 X 10 5 Tg" 2/3 exp(-6.12/Tg 1/3 ) 


2.56 


3.006 


4 Hex( 3 He, 7 ) 7 Bex 


3.2 X 10 6 T g " 2/3 exp(-9.72/Tg 1/3 ) 


2.56 


2.567 


4 HpW 6 T,i <v") 10 B v 
nl ~X \ !) D X 


2 4 x io 6 T _2/ ' 3 Mnf- Ifi 05/T 1 ^ 3 "! 


6.22 


6.154 


9 Bex(p, 6 Li) 4 Hcx 


1.9 x 10 11 Tg" 2/3 exp(-8.83/Tg 1/3 ) 


0.181 


0.991 


4 Hex(a,7) 8 Be x 


2.9 x 10 6 T g " 2/3 exp(-10.70/Tg 1/3 ) b 


7.89 


0.966 


6 Lix(o, 7 ) 10 Bx 


3.0 x 10 6 Tg" 2/3 exp(-16.98/T 9 1/3 ) 


3.38 


5.692 


7 Lix(Q, 7 ) n Bx 


2.7 x 10 7 Tg" 2/3 exp(-16.98/T g 1/3 ) 


7.89 


9.882 


7 Bex(a,7) n Cx 


6.6 X 10 7 T 9 " 2/3 exp(-22.26/T 9 1/3 ) 


7.89 


8.934 


4 Hex (d,X~) 6 Li 


2.37 x 10 8 (1 - 0.34T 9 )T 9 " 2/3 cxp(-5.33/T g 1/3 ) 


0.192 


1.131 


6 Bex(,e+^e) 6 Lix 


T 1/2 = 2.3 s 




3.349 


8 Bx(,e+^) 8 Bex 


T 1/2 = 1.2 s 




17.030 



a For nucli des a = i,j,k,... w ith mass numbers A a and spins g a , the reverse coefficients are defined 
as follows (Fowler et al. 1967) on the assumption that X~ particle is much heavier than nuclides: 

0.9867(gigj/s fc )A 3/2 for a process ix(j,l) k x, (9i9j I '(9k9l))( A j I 'A k ) 3/2 for a process ix(j,k)lx, and 
1.0134( 9i9i /(g fe9i5m ))(A J /(A fe A i )) 3/2 for a process ix(j,kl)™x- 
b The iS-factor for the reaction 3 He(o,7) 7 Be was used 



